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Description 

[0001] TTie present invention relates to mutti-slice helical computerized tomography and more particularly to an al- 
gorithm, method and apparatus for using the same which reduces the data acquisition time and data processing time 

5 required to generate an image. 

[0002] In computerized tomography (CT) X-ray photon rays are directed through a patient toward a detector. Atten- 
uated rays are detected by the detector, the amount of attenuation indicative of the make up (e.g. bone, flesh, air 
pocket, etc.) of the patient through which the rays traversed. The attenuation data is then processed and back-projected 
according to a reconstruction algorithm to generate an image of the patient's internal anatomy. Generally, the "back 

10 projection" Is perfonned in software but, as the name implies, is akin to physically projecting rays from many different 
angles within an image plane through the image plane, the values of rays passing through the same image voxels 
being combined in some manner to have a combined effect on the voxel in the resulting image. IHereinafter the data 
corresponding to rays which are back projected will be referred to as back projection rays. 

[0003] During data acquisition, if a patient moves, artifacts can occur in the resulting image which often render images 
15 useless or difficult to use for diagnostics purposes. For this and other reasons, as In other Imaging techniques, the CT 
Industry is constantly trying to Identify ways to reduce the duration of acquisition periods without reducing the quality 
of the data acquired. 

[0004] In addition, because huge amounts of data are acquired during an acquisition period and the processing 
methods for image reconstruction from the gathered data are relatively complex, a huge number of calculations are 

20 required to process data and reconstruct an image. Because of the huge number of required calculations, the time 
required to process collected data and reconstruct an image is appreciable. For this reason the CT industry is also 
constantly searching for new processing methods and algorithms which can speed up the reconstruction process. 
[0005] Various CT system features and procedures have been developed to increase data acquisition speed and to 
speed up the reconstruction process. Some of the more popular features and procedures including fan beam acquisi- 

25 tion, simultaneous multiple slice acquisition, helical scanning and half-scanning. In fan beam acquisition the source is 
collimated into a thin fan beam which is directed at a detector on a side opposite a patient. In this manner, a complete 
fan beam projection data set is instantaneously generated for a beam angle defined by a central ray of the source fan 
beam. The source and detector are rotated about an image plane to collect data from all (e.g., typically 360 degrees) 
beam angles. Thereafter the collected data is used to reconstruct an image in the image plane. Thus, fan beam ac- 

30 quisition reduces acquisition period duration. 

[0006] With respect to half-scanning, assuming a patient remains still during a data acquisition period, conjugate 
data acquisitions (i.e., data acquired along the same path from opposite directions) should be identical. In addition, 
using a fan beam, at least one ray can be directed through an image plane from every possible beam angle without 
having to perfomi a complete rotation about the patient. 

35 [0007] For example, referring to Fig. 3, an annular gantry opening 70 Is illustrated with a patient slice 42 disposed 
(support table not illustrated) therein and with respect to a Cartesian coordinate system where the Z-axis is into the 
Figure and defines a transport axis. A source 10 is illustrated in first, second, third and fourth positions as 90, 90\ 90" 
and 90'", respectively. When in the first position, source 10 generates a fan beam 40 which includes a central ray Rc 
and additional rays diverging therefrom along fan angles, the maximum fan angle being r. The beam angle B is defined 

40 as the angle formed by central ray Rc with respect to the vertical Y-axis. 

[0008] When in the fourth position, source 10 generates a fan beam 40"' which also includes a central ray (not 
illustrated) and rays diverging therefrom to fomi the fan beam. By rotating the source from the first to the fourth position 
in a clockwise direction data is collected at least once from every possible beam angle through slice 42 (i.e., the image 
plane). As known in the industry, data corresponding to every beam angle corresponding to a single image plane can 

45 be collected after a (n + 2r)/2n rotation about the patient. Because less than an entire rotation about the image plane 
is required to acquire the imaging data these acquisition methods and systems are generally referred to as half-scan 
methods and systems. Thus, half-scan acquisition has been employed to reduce acquisition period duration in con- 
junction with single row detectors. 

[0009] In addition, because relatively less data has to be processed in the case of half-scan imaging methods and 
so systems to generate an image, half-scan methods and systems also have the advantage of potentially reducing data 
processing and reconstruction times. 

[0010] While fan beams and half-scans have several advantages, often, during a diagnostics exercise a system user 
typically will not know the precise location within a patient of an object, cavity, etc. of Interest to be imaged. For this 
reason, it is advantageous for a system user to be able to generate several cross sectional images in rapid succession 
55 by selecting different image/reconstruction planes. In these cases rapid data processing is extremely important to 
minimize delays between image generation so that a user does not lose her train of thought between image views. 
[0011] Single slice detectors, fan beams and half-scans can be used to generate data in several different parallel 
image planes which, after data acquisition, can be used by a processor to generate an image anywhere between the 
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image planes through interpolation/extrapolation procedures known in the art. For example, assume that during two 
data acquisition periods first and second data sets were acquired which correspond to first and second parallel acqui- 
sition planes, respectively, the planes separated by 0.25 inches. If a user selects an image plane for reconstructing an 
image which resides between the first and second acquisition planes, interpolation between data In the first and second 

5 sets can be used to estimate values of data corresponding to the selected image plane. For Instance, assume that, 
among other rays, during the acquisition periods a first ray and a second ray were used to generate data in the first 
and second sets, respectively, and that the first and second rays were parallel (i.e. had the same beam and fan angles). 
In this case, by interpolating between the data acquired from the first and second rays generates an estimated value 
corresponding to a hypothetical back projection ray which is parallel to the first and second rays and which is within 

10 the image plane. By perfomiing such Interpolation to generate back projection rays for every beam and fan angle 
through the image plane a complete data set corresponding to the image plane is generated. 
[0012] White such systems work, unfortunately, the acquisition time required to generate data corresponding to many 
image planes is excessive and inevitable patient movement often causes image artifacts. 

[0013] One way to speed up data acquisition corresponding to several image planes is by employing a multi-row 

15 detector with a fan beam. In multi-row detector systems, a relatively thick fan beam is collimated and directed at a 
multi-row detector with a patient there between, each detector row in effect gathering data for a separate "slice" of the 
thick fan beam along the Z or translation axis perpendicular to a fan beam width. Despite each detector row having a 
thickness, in these systems it is assumed that the detected signals in each row correspond to a plane centered within 
the row as projected onto the isocenter Z. Hereinafter the central plane through a row will be referred to as a row center 

20 [0014] After data acquisition an interface enables a system user to select an image plane from within the area cor- 
responding to the collected data. The selected image plane is between the row centers of at least two adjacent detector 
rows. After image plane selection, a processor interpolates between data corresponding to adjacent rows to generate 
back projection rays corresponding to the selected image plane. When another image corresponding to a different 
Image plane is desired, after selecting the plane, the processor again identifies an acquired data subset for interpolation, 

2s additional processing and back projection. Thus, multi^row detector systems further reduce data acquisition period 
duration where several image planes may be selected for reconstruction. 

[0015] One limitation with multi-row detectors is that, during a single acquisition period, data can only be collected 
whteh corresponds to the detector thickness. To collect additional data corresponding to a greater patient volume, after 
one acquisition period corresponding to a first volume, the patient has to be moved along a translation axis until a 

30 second volume which is adjacent the first volume is between the source and detector. Thereafter a second acquisition 
process has to be performed. Similarly, to collect additional data corresponding to a third volume the patient has to be 
transported to another relative location with respect to the source and detector. Required translation without acquisition 
necessarily prolong the acquisition period and the additional acquisition time and aligning processes inevitably result 
in relative discomfort, additional patient movements and undesirable image artifacts. 

35 [0016] Helical scanning systems have been developed so that data can be collected during a single acquisition period 
without halting patient translation during the acquisition period. In a helical scanning system, the source and detector 
array are mounted on opposing surfaces of an annular gantry and are rotated there around as a patient is transported 
at constant speed through the gantry. The X-ray beam sweeps a helical path through the patient, hence the nomen- 
clature "helical scanning system". Data acquisition can be sped up by increasing operating pitch (i.e.. table translation 

40 speed relative to gantry rotation rate). After data Is acquired the data is processed to generate back projection ray 
estimates and account for data nuances which are caused by the helical acquisition. 

[0017] Various combinations of the fan-beam, multi-slice, half-scan and helical scanning features have been com- 
bined to realize synergies and have been somewhat successful. For example, one system combines a multi-row fan 
beam detector and a fan beam source with a helical scanning procedure to rapidly acquire Imaging data using a high 

^ pitch/high speed mode. For example, an exemplary system including a four row detector may support a 6:1 pitch 
wherein the detector crosses a reconstruction plane in 0.67 rotations. It is envisioned that an eight slice scanner will 
be able to support an 11 :1 pitch so that the detector crosses a reconstruction plane in 0.73 rotations. 
[0018] Referring again to Fig. 3, in these high speed helical scanning systems, during acquisition data is acquired 
with source 10 at position 90, the source and detector are rotated (while data is collected) about gantry opening 70 as 

so the patient 42 is transported there through. A processor collects data during transport and rotation from many different 
beam and fan angles. After source 10 rotates through a complete rotation and reaches position 90 again, additional 
data is gathered at that position. Because the patient 42 is transported along the Z axis during acquisition, while source 
1 0 is at the same location 90 relative to opening 70 at the beginning and at the end of the rotation, the source and data 
collected are at a different Z location relative to patient 42. Hereinafter data collected for the same beam and fan angles 

S5 but at different Z locations will be referred to as consecutively collected data. 

[0019] Cunrent interpolation techniques interpolate between consecutively collected data (i.e., data from source 10 
at the same beam angle (e.g., position 90 in Fig. 3) and fan angle but at different Z (i.e., translation axis) locations. In 
other words, current interpolation techniques require data from more than a single source rotation to generate an image. 
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In addition, because data from more than one rotation is required to interpolate, the data collection is relatively large 
and processing and reconstruction period durations are excessive. Moreover, because interpolation is between con- 
secutively collected data, the resulting image has a "thickness" characteristic which corresponds to a relatively thick 
patient volume which Is unsuitable or at least not optimal for many diagnostic purposes. 
5 [0020] An embodiment of the Invention will now be described, by way of example, with reference to the accompanying 
drawings, in whbh: 

Fig. 1 Is a perspective view of a CT apparatus used to practice the present Invention which Includes a detector 
array having rows and columns of detector elements and fan beam source; 

10 

Fig. 2 is a block diagram of CT control system which may be used to control the CT apparatus of Fig. 1 and which 
is useful for the purposes of practicing the present invention; 

Fig. 3 is a schematic view illustrating a gantry opening with a radiation source positioned at different locations with 
15 respect to the opening and also illustrating fan beams which con-espond to the beginning and end source positions 

of an exemplary half-scan; 

Fig. 4 is a schematic edge view illustrating a four row detector and a single image plane, the image plane illustrated 
in four distinct positions relative to the detector; and 

20 

Fig. 5 is a flow chart illustrating an exemplary embodiment of the present invention. 
A. HARDWARE 

25 [0021] Referring now to Fig. 1 , a CT scanner for use with the present invention includes a gantry 20 having an 
opening (70 in Fig. 3) supporting an x-ray source 10 oriented to project a fan beam 40 of x-rays along the beam axis 
41 through a patient 42 to a supported and opposed detector an-ay 44. The gantry 20 rotates to swing the beam axis 
within a gantry plane 38 defining the x-y plane of a Cartesian coordinate system. Rotation of the gantry 20 is measured 
by beam angle p from an arbitrary reference position within the gantry plane 38. 

30 [0022] A patient 42 resets on a table 46 which may be moved along a translation axis 48 aligned with the Z-axis of 
the Cartesian coordinate system. Table 46 crosses gantry plane 38 and is radiotranslucent so as not to interfere with 
the Imaging process. 

[0023] The x-rays of the fan beam 40 diverge from the beam axis 41 within the gantry plane 38 across a transverse 
axis 50 generally orthogonal to both the beam axis 41 and the translation axis 48 at a fan beam angle (. The x-rays of 
35 beam 40 also diverge slightly from the beam axis 41 and the gantry plane 38 across the translation axis 48. Referring 
also to Fig. 3, a maximum beam angle y is identified by symbol V. 

[0024] After passing through patient 42, the x-rays of the fan beam 40 are received by detector array 44 which has 
multiple columns of detector elements 1 8'. The detector elements 1 8' are an-anged in rows extending along the traverse 
axis 50 and columns extending along the translation axis 48. The surface of detector array 44 may be planar or may 

40 follow a section of a sphere or cylinder having a center at focal spot 26 or altematively at the system isocenter 

[0025] The detector elements 18' each receive x-rays and provide Intensity measurements along separate rays of 
the fan beam 40. Each intensity measurement describes the attenuation via a line integral of one fan beam ray passing 
through a portion of volume 43 of patient 42. In a preferred embodiment, volume 43 is greater than the slice volume 
measured by a conventional single slice fan beam CT system and the width of the detector array 44 is measured along 

45 Its columns. 

[0026] In a preferred embodiment, volume 43 is typically larger than the slice volume measured by a conventional 
fan beam CT system, and the width of the detector an-ay 44 as measured along Its columns is typically larger than the 
width of a single slice detector. The rows of detector elements 1 8' subdivide the fan beam detector array along the Z-axis. 
[0027] Referring now to Fig. 2, an exemplary control system for controlling the CT Imaging system of Fig. 1 has 

so gantry associated control modules 52 which include an x-ray control 54, a gantry motor control 56, a data acquisition 
system 62 and an imagery constructor 68. The x-ray control 54 provides power and timing signals to the x-ray source 
1 0 to turn it on and off as required under the control of a computer 60. The gantry motor control 56 controls the rotational 
speed and position of the gantry 20 and provides Infonnation to the computer 60 regarding gantry position. The data 
acquisition system 62 samples and digitizes intensity signals from the detector elements 18' of detector array 44 and 

55 the imagery constructor 68 receives the sampled and digitized intensity signals from the data acquisition system 62 
each identified as to column and row of the detector element of the detector array 44, and combines the intensity 
signals from the detector elements 1 8' according to the present Invention, and perfomis high speed Imagery construc- 
tion according to methods known in the art. 
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[0028] Each of the above modules is connected to Its associated elennents on the gantry 20 via slip rings 64 and 
serves to interface processor or computer 60 for performing various gantry functions. Slip rings 64 permit gantry 20 to 
rotate continuously through angles greater than 360° to acquire projection data. 

[0029] The speed and position of table 46 along the translation axis 48 Is communicated to and controlled by computer 
5 60 by means of table motor control 58. In addition, computer/processor 60 runs a pulse sequencing program to perform 
the inventive data processing method as described in more detail below. The computer 60 receives commands and 
scanning parameters via operator console 65 which Is generally a CRT display and keyboard. Console 65 allows an 
operator to enter parameters for controlling a data acquiring scan and to display the reconstructed image and other 
infomation from the computer 60. A mass storage device or memory 66 provides a means for storing operating pro- 
10 grams for the CT imaging system, as well as image data for future reference by the operator. Both computer 60 and 
the image constructor have associated electronic memory (not shown) for storing data. 

[0030] In operation, the gantry motor control 56 brings the gantry 20 up to a rotational speed and the table motor 
control begins translation of the table 46. The x-ray control 54 turns on the x-ray source 10 and projection data is 
acquired on the continuous basis. At each beam angle the projection data acquired comprises intensity signals 
15 corresponding to each detector element 1 8' at each particular column and row of array 44. 

B. THEORY 

[0031] It has been recognized that, In a multi-row data acquisition system which is run at high speed (e.g., see Figs. 

20 1 and 2), Instead of interpolating between consecutively collected data to estimate a back projection ray value within 
an image plane, interpolation can be perfonned between multi-row data generated with a source and a detector at a 
single location to achieve a number of advantages including reduction in acquisition period duration and reduction in 
data processing time after acquisition. Hereinafter, data acquired by different rows of a detector with the detector and 
source at a single position with respect to a gantry will be referred to as simultaneously collected data. 

25 [0032] With respect to reducing acquisition time, interpolation between simultaneously collected data enables image 
reconstruction using data corresponding to a "half-scan" (e.g., less than an entire source rotation). Herein the term 
"half-scan" is used to refer to any scan including less than a full source rotation and may, for example, Include 0.60 
rotations, 0.80 rotations, 0.74 rotations and so on. A half-scan rotation requires only a fraction of the time required to 
perfonn a complete rotation and therefore acquisition time is appreciably reduced. Because acquisition period duration 

30 is reduced, likelihood of patient movement and hence image artifacts is also appreciably reduced. 

[0033] With respect to reducing data processing time, the processing time is at least in part dependent upon the 
amount of data which must be processed to generate an image. As indicated above, with the present invention an 
image can be reconstructed using data corresponding to a fraction of a full rotation of the source about the patient as 
opposed to previous high speed helical systems which required data corresponding to more than a full source rotation 

35 about the patient. Clearly image processing time is appreciably reduced (e.g. , processing time is almost halved in some 
cases). By reducing processing time a more user friendly Interface which processes and displays Images more quickly 
can be configured. 

[0034] It has also been recognized that, in some systems which support only certain pitches or translation speeds 
relative to rotation speeds, the data generated during a high speed half-scan is insufficient for the purposes of image 

40 processing via Interpolation of simultaneously collected data (I.e., from data collected with the source in a single location 
relative to the patient. To this end, the data set con^esponding to a specific Image plane and a specific beam angle Is 
said to be a "complete data set" when a detector array (e.g., 44 in Fig. 1) simultaneously collects data on either side 
of the image plane. Similarly, the data set corresponding to a specific image plane and beam angle is said to be an 
"incomplete data set" when a detector array does not simultaneously collected data on either side of the image plane 

45 (i.e., for the plane and angle the anray only collects data on one side of the image plane). 

[0035] For example, assume that given a maximum fan beam angle r, data must be collected for at least 0.65 source 
rotation about a patient for half-scan image reconstruction, the 0.65 rotation corresponding to 0.325 rotation on either 
side of the image plane. Also assume that at a high speed pitch supported by the system, the width of a four row 
detector (where the width is between edge detector row centers) crosses the plane of reconstruction in 0.60 source 

so rotations about a patient. In this case, for the 0.3125 through 0.325 rotation before and the 0.3125 through 0.325 
rotation after the image plane, the acquisition will provide incomplete data sets (i.e., will not provide simultaneously 
collected data corresponding to points on both sides of the image plane) and therefore interpolation between simulta- 
neously collected data cannot be performed. Instead, for the 0.31 25 through 0.325 rotations before and after the Image 
plane, simultaneously collected data Is extrapolated to estimate back projection values. 

55 [0036] Referring now to Fig. 4, a schematic diagram illustrating an exemplary four row detector 44 is illustrated 
wherein the rows extend into the Figure. The four rows are identified by numerals 72, 74, 76 and 78, respectively, and 
together define a thickness between detector edges 100 and 102. Detectors 72 and 78 are edge detectors because 
they are on the ends of detector 44. Each row 72, 74, 76 and 78 has a row center 80, 82, 84 and 86, respectively 
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[0037] When data Is collected by a detector row, processor 60 (see Fig. 2) assumes that the data Is all collected at 
the corresponding row center as opposed to throughout the thickness of the row. For example, with respect to row 72, 
processor 60 assumes all data was collected at row center 80 and for row 74 processor 60 assumes all data was 
collected at row center 82 and so on. Referring to Figs. 1 and 4, detector 44 Is positioned on gantry 20 such that rows 

5 72 through 78 are orthogonal to translation axis 48 and extend along transverse axis 50. 

[0038] In some cases, prior to data acquisition, a system user will have performed some perfunctory exploratory 
imaging and hence will generally know the location of a desired image plane with respect to a patient* s anatomy. In 
this case the user could limit data acquisition to a data set which includes just enough data (e.g., data corresponding 
to a half-scan) for reconstructing an image corresponding to the desired plane. In other cases the user may not know 

10 the exact image plane desired and therefore may want to acquire data corresponding to a patient volume (e.g., 43 in 
Fig. 1) from which subsets of data can be extracted for reconstructing images in specific image planes. 
[0039] At this point it will be assumed that prior to acquiring imaging data, the actual plane through volume 43 which 
will be selected for imaging purposes is unknown. In this example, the image plane through volume 43 is selected 
subsequent to data acquisition via console 65 after which a specific set of the acquired data is selected for image 

15 processing. Nevertheless, for the purpose of simplifying this explanation, it will be assumed that we know the location 
of aspectfic image plane during data acquisition. This assumption enables oneto easily visualize the spatial relationship 
between the specific image plane and relative to detectors 72 through 78. 

[0040] Referring still to Figs. 1 , 3 and 4, during data acquisition, as patient 42 is transported through gantry opening 
70 a patient slice corresponding to the image plane moves along the Z axis and relative to detector 44. In Fig. 4 an 

20 exemplary image plane Pi is illustrated in four separate positions with respect to detector 44, the four separate positions 
identified by numerals 90, 90', 90" and 90"'. The four positions 90, 90'. etc., con'elate to four separate snapshots In 
time during a single data acquisition period as detector 44 is rotated about gantry opening 70 (see also Fig. 3) and 
patient 42, including the patient slice corresponding to image plane Pi. is translated through opening 70. In this example. 
Image plane positions 90, 90', 90" and 90"' in Fig, 4 con-espond to like numbered source positions in Fig. 3. 

25 [0041] Referring still to Figs. 3 and 4, clearly during data acquisition and as the image plane Pi traverses across the 
detector thickness (I.e. between edges 100 and 102), during some instances and corresponding beam angles, image 
plane Pi Is between adjacent detector row centers and during other instances and corresponding beam angles, image 
plane Pi is not between adjacent detector row centers. For example, on one hand, when source 10 and image plane 
Pi are at position 90, plane Pi is not between two detector row centers. On the other hand, at position 90' plane Pi is 

30 between row centers 84 and 86. Similarly, at position 90" plane PI is between row centers 80 and 82 and at position 
90"' plane Pi is not between adjacent row centers. Thus, for the data corresponding to positions 90' and 90", there is 
sufficient data for processor 60 to interpolate between adjacent detector row data to estimate back projection rays for 
the image plane Pi and source positions 90' and 90" (see Fig. 3) while at positions 90 and 90"' interpolation between 
simultaneously collected data is not possible and extrapolation is required. 

3s [0042] Referring to Fig. 2, after Imaging data required to generate an image within an image plane PI is acquired 
and stored, to generate an image, processor 60 performs the interpolation and/or extrapolation processes on a half- 
scan data set corresponding to plane Pi generating the estimated back projection data as described above. Thereafter 
half-scan weights are applied to the estimated back projection data to generate weighted back projection data and the 
weighted data is then combined via any of several well known back projection techniques to generate an image. The 

40 half -scan weights can be derived using several different half -scan algorithms. To this end, refer to the article entitled 
"Optimal Short Scan Convolution for Fanbeam CT' which was authored by D.L. Parker and which was published in 
Med. Physics Volume 9(2) in March, 1982. Nevertheless, the invention contemplates at least one prefen-ed algorithm 
wherein the half-scan weights Whs ^^r scan data acquired with the source angle between 0 and 2k and centered on 
n are given by solving the following equation: 



45 
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In this case a weight smoothing transfomnation may also be applied to the data according to the following equation: 

f(x)=3x^-2x^ 

withx = WHs(P.Y). (2) 

[0043] Alternative half-scan weight algorithms are contemplated. 

10 C. INVENTIVE METHOD 

[0044] Referring now to Fig. 5, an exemplary and preferred method according to the present Invention is illustrated. 
Referring also to Fig. 2, at process block 150, it is assumed that processor 60 has been preprogrammed with two 
separate numbers Including the number (#rhsr) rotations required for half-scan reconstruction for the specific system 
*5 configuration used with processor 60 and the number {#^^) of rotations required for detector 44 to cross a single image 
plane given a specific pitch supported by the system. With respect to the number of rotations required for haif-scan 
reconstruction tf^^sp that number is detenmined by solving the fottowing equation: 



20 



#,hs, = {B+2r)/(2B) (3) 



where r is the maximum fan beam angle (see Fig. 3). 

[0045] Continuing, at process block 152, processor 60 controls the system of Figs. 1 and 2 to acquire high speed/ 
high pitch data for several helical rotations and corresponding to volume 43 of patient 42. The acquired data Is stored 
in storage unit 66. At process block 1 54, a system operator uses console 65 to indicate one of a plurality of different 
image planes through volume 43 for which an image should be generated. 

[0046] Next, at block 156 processor 60 accesses data con-esponding to #rhs/2 rotations (I.e., corresponding to half 
the data required for half-scan reconstruction of either side of the Image plane PI) on either side of the Image plane Pi 
from mass storage unit 66. At decision block 158, processor 60 determines whether or not the number of Uj^^^ of 
rotations required for half -scan reconstruction is greater than the number U^^^ of rotations required for the detector 44 
to cross the image plane. If the number #^3,. of rotations required for half-scan reconstruction is not greater than the 
number #rcr of rotations required for the detector to cross the image plane PI, control passes to block 160. For example, 
if the number of rotations required for half-scan reconstruction Is 0.65 rotations and the number of rotations required 
for detector 44 to cross the image plane PI Is 0.73, all data required for generating an Image in the specified image 
plane can be determined via Interpolation and control passes to block 160. At block 160, processor 60 interpolates 
between simultaneously collected data for all data corresponding to rotations on either side of PI and for each 
fan beam angle p and ray angle y. After interpolation at block 160, control passes to block 162. 
[0047] Referring again to decision block 158, if the number of rotations required for half-scan reconstmction Is 
greater than the number of rotations required for the detector to cross the image plane Pi, control passes to block 
^ 1 68. In this case, interpolation can be used to identify data for generating an Image in an image plane corresponding 
to certain beam angles while extrapolation must be used to generate data for fonnlng the image from other beam 
angles. To this end, at process block 168 processor 60 Interpolates between simultaneously collected data for all data 
corresponding to #^/2 rotations on eltherside of image plane Pi. Then, at block 1 70 processor 60 extrapolates between 
simultaneously collected data for all data within #rhs/2 rotations on either side of the image plane Pi and outside #^0/ 
2 rotations on either side of the image plane Pi. In other words, extrapolation Is used to estimate projection values 
within the image plane only where interpolation Is Impossible. After extrapolation at block 1 70, control passes to block 
162. 

[0048] After the interpolation step at block 160 or the interpolation and extrapolation steps at blocks 168 and 170, 
at block 162 processor 1 60 calculates half-scan weights for each estimated backprojection value generating weighted 
values by solving equations 1 and 2 above. In the alternative, other half-scan weight algorithms are contemplated. 
After the weighted values have been generated, at block 164 processor 60 backprojects the weighted value data to 
generate an image and display the image on operator console 65. Until another command is received from the system 
operator at block 166, control continues to loop back up to block 164 and the image is displayed. At block 166, when 
a system operator uses console 65 to select another image plane Pi, control passes to block 154 and the process 
loops once again through the lower portion of Fig. 5 until an Image con^esponding to the newly selected image plane 
Pi is generated and displayed at block 1 64. 

[0049] It should be understood that the methods and apparatuses described above are only exemplary and do not 
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limit the scope of the invention, and that various nnodlfications could be nnade by those skilled In the art that would fall 
under the scope of the invention. For example, as indicated above, in some cases a system will know the precise 
location of a desired image plane with respect to a patient's anatomy. In that case only the data required to generate 
an image within the desired Image than would have to be acquired and the steps of selecting an image plane (i.e., 154 
5 In Fig. 5) and selecting specific data (i.e., 156 in Fig. 5) for reconstnjctlon would be avoided. In addition, referring again 
to Fig. 5, while the Inventive method Is described as Including comparison of rotation numbers at block 158, other 
methods for detenmlning If extrapolation will be used in addition to interpolation are contemplated. For example, systems 
which make this detenmination based on Z-locatlon are contemplated. 

[0050] For completeness, various aspects of the invention are set out in the following numbered clauses: 

10 

1. A method to be used with a CT system including a fan beam source (26) and a multi-row detector (44), the 
source (26) and detector (44) arranged on opposite sides of a transport axis (48), the detector having first and 
second edges (100, 102) separated along the transport axis and having separate detector rows (72, 74, 76, 78) 
parallel to the edges, the source directing a fan beam (40) including a central ray (R^) at the detector, each detector 

IS row subtended by a separate beam slice, the slices diverging from the central ray and including first and second 

edge slices which define beam edges orthogonal to the transport axis and which diverge from the central ray along 
a maximum fan angle r, the method for generating a cross sectional image of an object along an Image plane and 
comprising the steps of, with the source generating the fan beam: 

20 rotating (152) the source and detector about the transport axis at a rotation rate; 

transporting (152) the object relative to the fan beam (40) at a transport rate relative to the rotation rate such 
that the Image plane traverses the first and second detector edges after a partial rotation of the source about 
the transport axis; 

25 

while rotating and transporting, collecting (152) data corresponding to at least (n + 2r)/4ii rotation and less 
than one half rotation on either side of the image plane; and 

processing (160, 168, 170, 172, 1 64) the collected data to generate the cross sectional image. 

30 

2. The method of clause 1 wherein the central ray (RJ sweeps a beam angle through the object during data 
collection and the collected data Includes complete data sets for at least a subset of beam angles, each complete 
data set including detector row data on both sides of the Image plane and which corresponds to adjacent detector 
rows (72, 74, 76, 78) and wherein the step of processing Includes, for each complete data set, interpolating (160) 

35 between data in the set to generate an estimated value of a ray along the beam angle in the image plane. 

3. The method of clause 2 wherein the step of Interpolating (160) includes linearly interpolating between the data 
corresponding to at least one detector row on either side of the image plane. 

40 4. The method of clause 2 wherein the step of Interpolating (1 60) Includes helically weighting the data corresponding 

to at least one detector row on either side of the image plane. 

5. The method of clause 2 wherein the step of interpolating (1 60) includes applying a sine helical weighting to the 
data con-espondlng to all of the detector rows on each side of the image plane. 

6. The method of clause 2 wherein the step of transporting (152) includes transporting at a rate such that the 
fraction is at least {n + 2T)l2n) and wherein the collected data includes a complete data set for every beam angle 
within the (n -I- 2r)/2n) rotation. 

so 7. The method of clause 2 wherein the step of processing further includes the step of, after Interpolating (152), 

applying (162) half-scan weights to the estimated values to generate weighted values and back-projecting (164) 
the weighted values to generate the image. 

8. The method of clause 7 wherein the half-scan weights (162) are detennined by solving the following equation: 

S5 

^Hs (P.v) = ^^^ir- Pw=(«^)-r # P # p=(Jc/2)+r - 2y 
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10 



55 



Whs(P. Y) =10 P" # P # = (371/2) - r- 2y 
^HS (P^T) = ^^1^^ P"" # P" # Psup = (371/2) + r 
and then applying the following weight smoothing transformation: 



f(x)=3x^-2x^ 



where x=Whs{P. Y)- 



9. The method of clause 2 wherein the fraction is less than (tc + 2r)/2Tc such that the collected data includes at 
least some incomplete data sets, an Incomplete data set including detector row data on only one side of the Image 

15 plane for a specific beam angle and wherein the step of processing further Includes, for each incomplete data set, 

extrapolating (1 70) using data In the set to generate an estimated value of a ray along the beam angle in the image 
plane. 

1 0. The method of clause 9 wherein the step of processing further includes the step of, after interpolating, applying 
20 (1 62) half-scan weights to the estimated values to generate weighted values and back-projecting (1 64) the weight- 
ed values to generate the image. 

11 . The method of clause 10 wherein the half-scan weights are detemnined by solving the following equation: 

w^s (P.Y) = ^^^^ Pw=(7^)-r # p # p-=(7i/2)+r - 2y 

Whs(P. Y) = 1 0 p" # p # P^ = (371/2) - r- 2y 
,0 ^HS (P.Y) = ^^g^ P"" # P' # Psup = (3^/2) + r 

and then applying the following weight smoothing transfonnation: 

35 f(x)=3x^-2x^ 
where x=Whs(P, Y)- 

12. An apparatus for use with a CT system including a fan beam source and a multi-row detector (44), the source 
40 (26) and detector (44) arranged on opposite sides of a transport axis (48), the detector having first and second 

edges (100, 102) separated along the transport axis and having separate detector rows (72, 74, 76, 78) parallel 
to the edges, the source directing a fan beam (40) including a central ray (RJ at the detector, each detector row 
subtended by a separate beam slice, the slices diverging from the central ray and Including first and second edge 
slices which define beam edges orthogonal to the transport axis and which diverge from the central ray along a 
45 maximum fan angle r, the apparatus for generating a cross sectional image of an object along an image plane 

and comprising, with the source generating the fan beam: 

a processor (60) running a pulse sequencing program to perfomri the steps of: 

50 rotating (152) the source and detector about the transport axis (48) at a rotation rate; 

transporting (152) the object relative to the fan beam (40) at a transport rate relative to the rotation rate such 
that the image plane traverses the first and second detector edges (100, 102) after a partial rotation of the 
source about the transport axis (48); 



while rotating and transporting, collecting (152) data corresponding to at least (n + 2r)/47c rotation and less 
than one half rotation on either side of the image plane; and 
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processing (160, 168, 170, 172, 164) the collected data to generate the cross sectional image. 

13. The apparatus of clause 12 wherein the central ray (R^.) sweeps a beam angle through the object during data 
collection and the collected data includes complete data sets for at least a subset of beam angles, each complete 
5 data set including detector row (72, 74, 76, 78) data on both sides of the Image plane and which corresponds to 

adjacent detector rows and wherein the processor (60) runs the program to perfonn the step of processing by, for 
each complete data set, interpolating (1 52) between data in the set to generate an estimated value of a ray along 
the beam angle in the image plane. 

10 14. The apparatus of clause 13 wherein the processor (60) runs the program to perfonn the step of Interpolating 

(152) by linearly interpolating between the data con-esponding to at least one detector row on either side of the 
image plane. 

15. The apparatus of clause 13 wherein the processor (60) runs the program to perform the step of interpolating 
15 (1 52) by helically weighting the data corresponding to at least one detector row on either side of the image plane. 

16. The apparatus of clause 13 wherein the processor (60) runs the program to perfomn the step of interpolating 
(152) by applying a modified sine helical weighting to the data corresponding to all of the detector rows on each 
side of the image plane. 

20 

17. The apparatus of clause 13 wherein the processor (60) runs the program to perform the step of transporting 
by transporting (1 52) at a rate such that the fraction Is at least {k + 2r)/2n) and wherein the collected data includes 
a complete data set for every beam angle within the {k + 2r)J2n) rotation. 

25 18. The apparatus of clause 13 wherein the processor (60) further mns the program to perfonm the step of, after 

interpolating (152). applying (164) half-scan weights to the estimated values to generate weighted values and 
back-projecting the weighted values to generate the Image. 

1 9. The apparatus of clause 1 8 wherein the half-scan weights (1 62) are determined by solving the following equa- 
30 tion: 

^HS (P'Y) = ^^p^ Pinf=(«/2)-r # P # P =(n/2)4.r - 2y 
35 Whs(P.y) = 1 .0 p- # p # P* - (3n/2) - r - 2y 

^Hsi^^y) = ^^^T^ P"" # P' # Psup = i^^f^) + r 

and then applying the following weight smoothing transfomnation: 

40 

f(x)=3x^-2x^ 

where x=Whs(P. Y). 

45 

20. The apparatus of clause 13 wherein the fraction is less than (n 4- 2r)/27i: such that the collected data includes 
at least some Incomplete data sets, an incomplete data set including detector row data on only one side of the 
image plane for a specific beam angle and wherein the processor (60) runs the program to perform the step of 
processing by, also, for each incomplete data set, extrapolating (1 70) using data In the set to generate an estimated 

50 value of a ray along the beam angle in the image plane. 

21 . The apparatus of clause 20 wherein the processor (60) runs the program to perform the step of processing by, 
also, after interpolating (1 52), applying (1 62) half-scan weights to the estimated values to generate weighted values 
and back-projecting the weighted values to generate the image. 

55 

22. The apparatus of clause 21 wherein the half-scan weights (1 62) are determined by solving the following equa- 
tion: 
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Pw=(n^)-r # p # P'=(ic/2)+r - 2y 




Whs(P.Y) = 10 



p* # P # p"" = (3n/2) - r - 2y 

Psup = (3'^/2) + r 



and then applying the following weight smoothing transfonnation: 



f(x)=3x^-2x^ 



where x=Whs(P. Y)- 



Claims 



1 . A method to be used with a CT system including a fan beam source (26) and a multi-row detector (44), the source 
(26) and detector (44) arranged on opposite sides of a transport axis (48), the detector having first and second 
edges (100, 102) separated along the transport axis and having separate detector rows (72, 74, 76. 78) parallel 
to the edges, the source directing a fan beam (40) Including a central ray (RJ at the detector, each detector row 
subtended by a separate beam slice, the slices diverging from the central ray and including first and second edge 
slices which define beam edges orthogonal to the transport axis and which diverge from the central ray along a 
maximum fan angle F, the method for generating a cross sectional image of an object along an image plane and 
comprising the steps of, with the source generating the fan beam: 

rotating (152) the source and detector about the transport axis at a rotation rate; 

transporting (152) the object relative to the fan beam (40) at a transport rate relative to the rotation rate such 
that the image plane traverses the first and second detector edges after a partial rotation of the source about 
the transport axis; 

while rotating and transporting, collecting (152) data corresponding to at least {n + 2r)/47i rotation and less 
than one half rotation on either side of the image plane; and 

processing (160, 168, 170, 172, 1 64) the collected data to generate the cross sectional image. 

2. The method of claim 1 wherein the central ray (R^.) sweeps a beam angle through the object during data collection 
and the collected data includes complete data sets for at least a subset of beam angles, each complete data set 
including detector row data on both sides of the Image plane and which corresponds to adjacent detector rows 
(72, 74, 76, 78) and wherein the step of processing Includes, for each complete data set, interpolating ( 1 60) between 
data in the set to generate an estimated value of a ray along the beam angle In the image plane. 

3. The method of claim 2 wherein the step of interpolating (160) includes linearly interpolating between, or helically 
weighting, the data con-esponding to at least one detector row on either side of the image plane. 

4. The method of claim 2 wherein the step of interpolating (1 60) includes applying a sine helical weighting to the data 
corresponding to all of the detector rows on each side of the image plane. 

5. The method of claim 2 wherein the step of transporting (1 52) includes transporting at a rate such that the fraction 
is at least {n + 2r)/2n) and wherein the collected data Includes a complete data set for every beam angle within 
the {n + 2T)f2n) rotation. 

6. An apparatus for use with a CT system including a fan beam source and a multi-row detector (44), the source (26) 
and detector (44) arranged on opposite sides of a transport axis (48), the detector having first and second edges 
(100, 102) separated along the transport axis and having separate detector rows (72, 74, 76, 78) parallel to the 
edges, the source directing a fan beam (40) including a central ray (R^) at the detector, each detector row subtended 
by a separate beam slice, the slices diverging from the central ray and including first and second edge slices which 
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define beam edges orthogonal to the transport axis and which diverge from the central ray along a maximum fan 
angle r, the apparatus for generating a cross sectional image of an object along an image plane and comprising, 
with the source generating the fan beam: 

5 a processor (60) running a pulse sequencing program to perfomi the steps of: 

rotating (152) the source and detector about the transport axis (48) at a rotation rate; 

transporting (152) the object relative to the fan beam (40) at a transport rate relative to the rotation rate such 
10 that the Image plane traverses the first and second detector edges (100, 102) after a partial rotation of the 

source about the transport axis (48); 

while rotating and transporting, collecting (152) data corresponding to at least (n + 2r)/4jt rotation and less 
than one half rotation on either side of the image plane; and 

15 

processing (160, 168, 170, 172, 164) the collected data to generate the cross sectional image. 

7. The apparatus of claim 6 wherein the central ray (RJ sweeps a beam angle through the object during data collection 
and the collected data includes complete data sets for at least a subset of beam angles, each complete data set 
20 including detector row (72, 74, 76, 78) data on both sides of the image plane and which corresponds to adjacent 

detector rows and wherein the processor (60) runs the program to perform the step of processing by, for each 
complete data set, interpolating (152) between data In the set to generate an estimated value of a ray along the 
beam angle in the image plane. 

25 8. The apparatus of claim 1 3 wherein the processor (60) runs the program to perform the step of interpolating (1 52) 
by linearly interpolating between, or helically weighting, the data corresponding to at least one detector row on 
either side of the image plane. 

9. The apparatus of claim 7 wherein the processor (60) njns the program to perfonm the step of interpolating (152) 
30 by applying a modified sine helical weighting to the data corresponding to all of the detector rows on each side of 

the image plane. 

10. The apparatus of claim 7 wherein the processor (60) runs the program to perform the step of transporting by 
transporting (152) at a rate such that the fraction is at least {n + 2r)/2n) and wherein the collected data includes 

35 a complete data set for every beam angle within the {n + 2r)/27c) rotation. 
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